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Abstract 

Homogenization appeared more than 100 years ago. It is an approach to 
study the macro-behavior of a medium by its micro-properties. In mathematics, 
homogenization theory considers the limitations of the sequences of the problems 
and its solutions when a parameter tends to zero. This parameter is regarded 
as the ratio of the characteristic size in the micro scale to that in the macro 
scale. So what is considered is a sequence of problems in a fixed domain while 
the characteristic size in micro scale tends to zero. But for the real situations 
in physics or engineering, the micro scale of a medium is fixed and can not be 
changed. In the process of homogenization, it is the size in macro scale which 
becomes larger and larger and tends to infinity. We observe that the homoge- 
nization in physics is not equivalent to the homogenization in mathematics up 
to some simple rescaling. With some direct error estimates, we explain in what 
means we can accept the homogenized problem as the limitation of the original 
real physical problems. As a byproduct, we present some results on the math- 
ematical homogenization of some problems with source term being only weakly 
compacted in H~ l , while in standard homogenization theory, the source term is 
assumed to be at least compacted in H . A real example is also given to show 
the validation of our observation and results. 
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1 Introduction 



Homogenization appeared more than 100 years ago. It is an approach to study the 
macro-behavior of a medium by its micro-properties. The origin of this word is related 
to the question of replacement of the heterogenous material by an "equivalent" homoge- 
nous one. The earliest papers dealing with the problem of this type are [1] , [2] , and the 
good survey of the results until 1925 is in [3]. The name of "homogenization" was first 
introduced by I. Babuska [I]. In physics, mechanics and engineering, homogenization 
is widely used to study the property of medium or material by the macro-behavior 
in stead of the complicated micro structure. The systematic mathematical theory of 
homogenization was built in [5]-[T5] and so on. But, does the mathematical theory 
describe the physics or engineering questions exactly? It seems hard to give a positive 
answer. In this paper we take the flow transport problem in the periodic heterogenous 
porous medium as an example to demonstrate the difference between the homoge- 
nization in mathematics and physics. Any other examples such as heat or electric 
conductivity and mass transfer will lead to the same conclusion. 

In mathematics, we consider the limitation of a sequence {w e } e >o such that 

-V • (A £ (x)Vu £ (x)) = f, in fi, 

(1.1) 

u £ = 0, on <9f2, 

where Q G R (d — 3) is occupied by the heterogenous porous medium and the per- 
meability coefficient A £ (x) = A(-) with A(y) being periodic with respect to y G Y = 
[0, l] d . u £ is the flow pressure in the medium and / is the source. Just for simplicity, 
we take Q = (0, l) d . By [3] and [T5], under the following assumption: 

(HO) : f G L 2 (VL) and there exist two positive constants A, A independent of e, 
such that AKI 2 < fK(y)£ < A|£| 2 , Vy, £ G R d , (1.2) 

we have that as e — > 0, there exists a u° G Hq(Q) such that u £ — ^ u° weakly in Hq(Q), 
where the u° is the solution of the homogenized problem 



1.3) 



-V • (A°Vu°(x)) = f, in Q, 

u° = o, on on, 

with A^j = rpr f Y {Aij + A ik ^-)dy and being the solution of the cell problem : 

-V v ■ (A(y)V y W(y)) = V y ■ (A(y) ej ), in Y = [0, l] d , 
N'(y) is Y-periodic and < N j >= ^ J Y N j (y)dy = 0. 1 ' ) 



Please note that in the above limit process, when the period size e tends to zero, 
the domain Q does not change, so there will be more and more e-periods contained in 



the whole domain. But this is not the case in the physics. In a physical or engineering 
problem, the size of the periodic micro cell structure can not be changed. What we 
consider in the physics is the following problem (take the flow in porous media as an 
example) 

-V • (K (f ) Vp{x)) = f, in n D = (0, D) d , 

(1.5) 

p = 0, on dflo, 

where l(m) and D(m) are the characteristic lengthes of the periodic micro cell and 
the whole medium respectively and I <C D, p (N/m 2 = Kg/m/s 2 ) is the pressure, 
K (m 3 s/kg) is the permeability coefficient and the source term / (1) = y with 
fo (kg/m 3 /s) and the density p (kg/m 3 ). 

We want to get the effective coefficient and the homogenized problem in physics 
and try to explain in what sense we can expect this result is valid. In the early work 
of I. Babuska [16] on the homogenization approach in engineering, he first pointed out 
that: '/ is a given parameter, with physical meaning which cannot be changed, e.g. 
cannot be made "sufficiently" small'. In his other work [TTj . |18j . he mentioned that e 
in fll.ip is the ratio of micro-scale (cell scale) to the macro-scale. If we set e = and 
I is fixed, then e — > means that D tends to infinity. The problem is if we take the 
transformation e = 4, can we transfer the physical problem (II. 5p to the mathematical 
problem (II. ip while the Assumption (HO) is still valid? 

Let us take a variable transformation x = mathematically, then the pressure p(x) 
satisfies that 

-V 2 • (D~ 2 K (§) Vgp(2c)) = fix), mfl = (0, l) d , 

(1.6) 

p{x) = 0. on dQ. 

It is obvious that the Assumption (HO) can not be satisfied even for a simple example 
as / = 1, since D~ 2 K(^) -»■ 0, as e -> 0. 

Let's go back to physics. What's the meaning of the variable transformation x — |j 
from (II. 5p to (II. 6p ? The only difference is the unit of length: in ( II. 5p . the length unit 
is one meter; in (11 .6p . the length unit is D meters. 

Therefore, if we take a variable transformation x = physically, then the pressure 
p{x) should satisfy the same Conservation Law of Mass, which is independent of the 
unit in length, 

-V • (K (I) Vflxj) =f, inQ = (0, l) d , 

(1.7) 

p — 0, on dVL, 

but all the physical quantities: the pressure p, the permeability coefficient K, the 
source term / and the gradient operator V are all measured in the new length scale, 



x = — (Dm), K = — ((Dm) 3 s / kg) , p = Dp (kg/Dm/s 2 ), 

V = DV(1/Dm), f = f(l/ s ). 

It should be noted that K = will tend to zero as D tending to infinity, since K 
is the fixed physical quantity in the original unit system. So the Assumption (HO) can 
not be valid either and we can not expect the validity of the homogenization theory 
even for a constant source term. 

So far, we see that the homogenization of the physical problem (jl.5p can not easily 
fall into the mathematical framework (ll.ip - (ll.4p by a direct transformation in mathe- 
matics or physics. In the following, we assume that D is sufficient large but fixed. We 
consider two different situations for the source term /. The first situation is that / has 
no micro- structure, and the other is that / has micro structure with period (0, l) d , for 
example, the source term may have the form as / = fi(x, x/l) + V ■ f'2(x, x/l). We first 
give the homogenized problem with the effective coefficient K° for different situations. 
We present the error estimate between the pressure p and the first order expansion p 1; 
then try to understand in what sense the homogenized problem is a limitation of the 
original problem. It is worthwhile to point out that so far the mathematical homoge- 
nization theory (for which we consider the limitation as I — Y 0) is still incomplete for 
the second situation, since the source term is only convergent weakly in if -1 (see [19] 
and [20] ). The homogenization theory for this kind of problem may have independent 
interests. 

The outline of paper is as follows: in §2 we discuss the situation that the source 
term has no micro- structure; in §3 we discuss the situation that the source term has 
micro-structure; in §4 a real example is given to show the validation of our observation 
and results. 

2 The source term has no micro-structure 

we say / has no micro-structure if / does not contain any micro-scale information at 
the scale comparable to or less than /. We discuss two different cases in this situation. 

• Case a: / G L 2 (R d ); Case b: / e L°°(i2 d ). 

By unit transformation, we get (11. 7p . which we have known do not satisfy Assump- 
tion (HO). In order to let the Assumption (HO) be valid, we introduce a new setting of 
problem as 

-V • (K £ (x)Vp £ (x)) = J(x), in tt = (0, l) d , 

(2.1) 

p £ = 0, on dfl, 



with 

' K e (x) = K(2) = D 3 K(I), ( K £ (x) = #(§ ) = D 3 K(^), 

a: ^ p £ (x) = D~2p(x), b: < p £ {x) = D~ 3 p(x), (2.2) 
7(2) = Dlf(x), { J(x) = fix). 

It is easy to check that Assumption (HO) is satisfied for the both cases. 

Then we know from homogenization theory (see [2~T] ) that there exists a po such 
that 

V s ~^ Po weakly in Hq(Q), 

(2.3) 

K £ Vp £ K°Vp weakly in {L 2 {Q)) d , 
where p is the solution of the homogenized problem of (12.11) : 

-V ■ (K^po(x)) =J, in = (0, l) d , 

(2.4) 

Po = 0, on dfl, 

with K9. = J Y (K ij (y) + K ik (y)^^)dy and iV J '(f ) = N j (y) solving of the cell problem: 



-V, ■ (K(y)V y (N^y) + Vj )) = 0, in Y = (0, l) d , 
N j is Y-periodic, J Y N j (y)dy = 0. 



(2.5) 



If we denote by 



p 1= p + £ N^, (2.6) 



then there exits a positive constant C independent of e such that 

\\P £ -Pi\\m m <Cei (2.7) 

After carefully checking the exact dependence of the constant C in ( 12. 7ft . we have 

Proposition 2.1 // i/ie coefficient and the source term of equation \2. 1\) satisfy the 
Assumption (HO), then there exists a positive constant C independent ofe, p £ , po, such 
that 

\\^{P £ -Pi)\\mn) < Ce^{\\G\\ L ~ {Y) + A\\N\\ L ~ {Y) )\\J{x)\\ L 2 (n) 

+C7e4(A)(i + ||iV|| i0 o (y) )||7(x)|| L2{ n), (2.8) 



A 2 

\\p £ -Po\\mn)< Ce±(\\G\\ L ~ {Y) + A|| N\\ L ^ Y) )\\J(x) \\ L2(U) . (2.9) 



Here G = (G 1 , ...,G d ) with G^ being a skew- symmetrical matrix satisfying (see \TBj) 
G{ k = K(y) tJ + K(y) ik — - K°, j = 1, d. (2.10) 



As to energy, we have 



Wp £ ■ K £ Wp £ - / Vp • K°Vp Q \ -+ 0, as e -»■ 0. (2.11) 
n Jn 

By the inverse transformation of (I2.2p and changing the unit in length from Dm to 
m, we obtain the following homogenized problem of (II. 5ft from ( 12.4ft and (12. 5 j) 



-V ■ (ir°Vp (x)) = /, in fi D = (0, D) 



(2.12) 



p = 0, on 8Qd, 
with K° determined as follows: 



Kl, = ^(^M + ^to)^)* 



<2 



where iV J ^yj = N^{y) (j = l...d) is the solution of cell problem 

-V, ■ {K{y)V y {W(y) + y 3 )) = 0, in Y = (0, 1 

is Y-periodic and J Y N^(y)dy = 0. 
It's easy to find that (I2.14p is equivalent to the following problem 
-V ■ (K(f )V(W(f ) + = 0, in Y = (0, l) d ., 

N j is Y t periodic and ^ Lnd N j (j)dx = 0. 



(2.14) 



(2.15) 



Remark 2.2 In fact, if we simply regard the micro size I in (jl.5p as a small parameter 
and formally apply the mathematical homogenization theory, we would obtain the same 
homogenization settings for ( 11.51) as (12. 121) - (12.151) . Furthermore we would still have the 
following mass balance or homogenization rule as : for any r\ e R d , if p v solves 

{ -V • (K(f ) Vp„) = 0, zn^ = [0,/] d , 

(2.16) 

«. Pr) — ^ ■ x is periodic in [0, 



then (Vp v ) = ja f Y Vp v dx = n and by A2.13\) and A2.15\) ,we have 

(K(j)V Pv ) Rd = (K(j)Vp v ) Yl = K°(V Pv ) Yl = K°(Vp v ) Rd , (2.17) 

which means the mass balance between the micro and macro scales. This is the reason 
why the homogenized coefficient K° is also called as the effective coefficient of K. 

The relationship between po(x) and po(x) are 

(Po(x) = D~^p (x), ( po(x) = D 2 p (x), 
b: I (2.18) 
K° = K°, [ K° = K°. 

By (12. 18)) and Proposition (12. ip . we can obtain the next theorem 

Theorem 2.3 If p is the solution of (II. 5ft . pi is defined as follows 

Pl =p + IN^. (2.19) 
and f G L 2 (R 3 ), then there exists a positive constant C independent of D such that 

i |V( ^S )|2 ^ " ^(^ll/^ll^ + ^ll/^ll^^)- (2-20) 

/ n J^-§| 2 ^<^) 2 H/(-)lli^)- (2-21) 
As to energy, there exists a positive const C independent of D, such that 

1-^3 / Vp-KVp-^ f Vp - K°V Po \ < cl, (2.22) 
u Jn D u Jn D u 

which means the convergence of the density of energy; 

Theorem 2.4 If p is the solution of ( 11.5)) . p\ is defined in (12.191) and f G L°°(R 3 ), 
then there exists a positive const C independent of D such that 



ji/ \ 2d ^ C ^nm\\Uu„y (2.24) 



D 



As to energy, we have 

1 I V-^ • KV^dx - — [ • K°V^dx\ -> 0, as D -» oo. (2.25) 

Jn D D D D Jn D D D 



1^3 



The above theorems explain in what sense we can accept the homogenized problem 



3 The source term has micro-structure 



In this situation, we will discuss the source term with the following micro-structure, 
-V • (K(f ) Vp(a;)) = / (x, f ) + V • F (x, f ) , in Q D = (0, D) d , 

(3.1) 

p — 0, on dViu- 
By unit transformation, we get 

-V • (K (I) Vp(x)) = / (x, § ) + V • F (x, f ) , in n = (0, 



(3.2) 



p — 0, on <9fi. 



Setting 



**® = (3.3) 



we obtain 

-V • (K e (x) Vj9 £ (x)) = / (x, f ) + V • F (x, f ) , in = (0, T 
p £ = 0, on dVl. 



d 

(3.4) 



The homogenization for this kind of problem may have independent interest, since the 
source term here is only weakly convergent in if -1 (f2) as e — > 0. The standard theory 
only treats the case that the source term is strongly convergent in if -1 (f2) see(|15]).In 
[T9] and [20] some incomplete results were present for this case. We will establish the 
homogenization theory for (13. 4p . 

In the beginning, we introduce an important lemma that will be used later. 

Lemma 3.1 If f(x,y) G L°°(Y, C 0,1 (fi)) and is Y-period with respect to y, where Q 
is an arbitrary bounded open subset of R d and Y = [0, l] d , then /(x, ~) — > Mf(x) = 
|F| Iy f( x iy)dy i n H~ 1 (Q) and there exists a constant C > such that 

\\f(x,-)-Mf(x)\\ B - Ha) <Ce\\f\\ L 

The proof of lemma is similar to the lemma 1.6 in |22j . 
Theorem 3.2 If u e is the solution of the following problem 

-V ■ (A £ (x)Vu £ ) = f (x, f ) + V ■ F (x, f ) , m SI, 

u £ = 0, on dQ, 



(3.6) 



where A e (x) = A (| ) is symmetric satisfying uniformly elliptic condition, i.e. there 
exist two positive constants X, A independent of e, such that A|£| 2 < C, T K(y)^ < 
A|£| 2 , Vy,£ G M. d , and A(y) is Y-period. f(x,y) and F(x,y) are bounded and Y-period 
with respect to y then as e — > 0, 



u e — ^ u° weakly in Hq(Q), 

A £ Vu £ + F(x, f ) A°Vu° + F°, weakly in (L 2 (tt)) d . 
u° is the solution of the following homogenized problem 

-V • (A°Vu°) = f°(x) + V • F°(x), in ft, 

u° = 0, on <9ft, 
with A , f°, and F° defined as follows 



(3.7) 



(3- 



A° j = (A ij (y)+A ik ^) Y , 

f°(x) = (f(x,y)) Y , (3.9) 
F°(x) = (F(x,y) + A(y)V y w(y)) Y , 



and Ni(y), w(x,y) solving the cell problems : 

-V y ■ (A(y)V y W(y)) = V y ■ (A(y)e j ), m Y, 

Ni(y) is periodic in Y, (iV J )y = 0, 

-V„ • (A(y)V y w(x,y)) = V y ■ (F(x,y)), in Y, 

w(x,y) is periodic in Y, (w)y = 0. 
Further more, if we denote by 



(3.10) 



(3.11) 



du° 

Ul = u ° + eN J — + ew, (3.12) 
dxj 

and assume thatA(y), f(x,y), F(x,y) are smooth enough and w(x,y) G lf 1,00 (fixF), 
then there exists a positive constant C independent of e such that 

||V(« e -Wi)|U 2(n) <Ce5 (3.13) 

Proof. By the standard asymptotic expansion method, we can get the equations (13.81) - 
fl3TT2|) ([19], [20]). We first use Tartar's method (|19J) to prove We denote by 

f £ = f(x, |) and F e = F(x, |) for short. By the regularity of elliptic equation, we 
obtain u e is bounded in Hq(Q) and £ e is bounded in (L 2 (Q)) d , where £ e = A e 'Vu £ + F e 
is a vector-function and satisfies 

/ i e Vvdx=(f,v) H -^ Hhmi Vve Hl(Q). (3.14) 
Jn 



By the compact property, there exists a subsequence (still denoted by e), such that 

u £ — ^ u° weakly in Hq(Q), 



e -» £°weakly in {L 2 {Q)) d , 
Taking e — > in (13.141) . we have 

/ ?Vvdx = (f,v) H - miB u n) , Vve H l (n). (3.16) 
Jn 

Therefor, (13.71) is proved if we show that 

£° = A°Vu° + F°. (3.17) 

If we set 

X £ x = \-x + eN x (^y (3.18) 

with A G R d and N\(y) solving the following problem 

-V, • (A(y)V y N x (y)) = V, • (A(y)X), in Y, 



N\ is periodic in Y and (N\) Y = 0, 

the we have the following limitation: 

X\ — A • x, weekly inif 1 (f2), 
X\ — > A ■ x, strongly inL 2 (f2). 



Introduce the vector function 



By the definition of x £ \ an d (I3.19p . we can easily obtain 

1)1 -± {A°f\ weakly in L 2 {{Q)) d , 
J n if x -Vvdx = 0, Vwe H%(Sl). 



(3.19) 



(3.20) 



r,l = A T V y xl- (3-21) 



(3.22) 



For any if G T>(Q), choose <px\ as the test function in ( 13. 14|) and <pu e as the test 
function in ( I3.22[) . We have 

/ f " Vxl^dx + [ e ■ Vcpxldx = (A ipxi) H -i(n), H i(n), V <p G V(Q), (3.23) 
Jn Jn 



I r)l ■ Vu e ipdx + / 7]\ ■ Vipu £ dx = 0, V if G V{VL). 
n Jn 



(3.24) 



By the definition of r] £ x , we have 

A £ Vu e ■ Vxl = (^ £ ) T Vxa ■ Vn £ = V { ■ Vu £ . 
From 0321-0321, we have 

/ e-V<pxidx+ [ F e -W X \vdx- [ V{-V<pu £ dx = (f £ ,<px x ) H -Hn),Hi(n)- (3-25) 
Jn Jn Jn 

Taking e -> 0, by and (|3T22]h we obtain 

I £° ■ Vy?(A ■ x)dx + / ^ / F(a;, y) • (A + V y N\)dy — [ X(A°) T ■ V^u°dx 
Jn Jn Jy Jn 

= (f°, (A • x)ip) H -i {n);H i {n) , (3.26) 
which can be rewritten in the form 

/ £° ■ V[<^(A • x))dx - I £° • V Xdx + I cpdx I F(x, y) ■ (A + V y N x )dy 
Jn Jn Jn Jy 

- / A(A°) T • Vym°dx = (/°, (A • x)<p) H - Hn)iHh{n) , (3.27) 
By (J3~H, we g et 

/ £°-^Ada; = f <pdx [ F(x,y) • (A + V y N x )dy + [ \{A°) T -Vu°cpdx. (3.28) 
Jn Jn Jy Jn 

If let A = ei, we can obtain 

£° = (A°Vu% + if. (3.29) 
Here we use the following relationship f[20j) 

i? = J(Fi(x,y) + (A(y)V v w)i)dy = J(F(x,y)(ei + V y N l ))dy. (3.30) 

Then we have proved (13. 7p . 

Next we will give the error estimate 

\\V(u £ - Ul )\\ L2(n) <Cel (3.31) 

Following the argument in [15] , we first compute 



(71 1 (Jl I 

A £ \7 Ul - A°\7u° = (A £ ( ej + \7 y N j ) - A° ej ) ^— + eA E N j V 



dxj dxj 



gKy) d ^- + eA £ N^{^ 

OXj OXj 



(3.32) 



with gi(y) = A £ (ej + V^iV- 7 ) — A°ej. By the cell problem ( I3.55P and the definition of 
A in ([33D, we find 

(9 j (y))Y = and V y -gi(y) = 0. (3.33) 
Then there exists a skew-symmetric matrix G- 7 ([IS]), such that 

g* = Vy&, Gi k = -Gi, and G\ k G H^ er {Y). (3.34) 

By this property, we can get: 

9 ^. £V .(^)-^V(^), (3.35, 
dxj \d Xj J \dxjj 

and 

- V • (A £ V(u £ - Ul )) = R £ + V-r £} (3.36) 
with R £ = V-(F(x, y) + A(y) V y w(x, y) - F°{x))+(f(x, y)-/ (x)), and r £ = eGl k £^+ 



sA £ ,N j 92u ° 



ij dxjdxk' 

If we introduce the boundary corrector 9 £ as 

-V ■ (A £ (x)V6 £ ) = 0, in fi, 



Q £ = -eN^-ew, on 



(3.37) 



then e £ = u £ — u\ — 9 £ G Hq(Q) satisfies 

- V- (A £ V(u £ - Ul -0 e )) = J R £ + V-r £ , (3.38) 

and the weak form is 

/ A £ Ve £ V(f)dx= I R £ cj)dx+ I r £ V<pdx, V0 G H^(Q) (3.39) 
in Jn Jn 

Since e e G Hq(CI), taking = e e , we get 

/ A £ Ve £ Ve £ dx= / R £ e £ dx+ / r £ Ve £ dx, V0 G ^(fi) (3.40) 
Jn Jn Jn 

By the elliptic condition and assumption that G G L°°(Y), the second term in the 
right hand can be estimated as follows 

d 2 u° 

||r £ ||i 2(n) <5(||G||^ (y) +A||iV||^ (y) )||^-|| L2(Q) . (3.41) 
By the definition of F°(x) in (13. 9p . we obtain 

V, • (F(x, y) + A{y)V y w{x, y) - F°{x)) = (v x ■ F - jf V x ■ F{x, y)dy^j 

+V X ■ (A(y)V y w(x, y)) ~ f / V* • (A(y)V y w(x, y))dy\ (3.42) 



It can be easily checked that 



R £ (x,y)dy = 0. 



(3.43) 



Y 



By lemma [37T1 and lemma 1.6 in [22], the first term in the right hand of (I3.40p can be 
estimated as follows: 

(3.44) 



/ R £ (x,y)e £ dx < Ce\\e e \\ H u n) 
Jn 



So we obtain 



1 d 2 u° 
|V(M e - Ml -r)|| i2(n) < Ce-(\\G\\ L ~ + A\\N\\ L ~)\\ ^ ^ Hfl^n) 

+Ce 



(3.45) 



Next we estimate || V# £ ||l2(q). Multiplying 6 s + e^fN 3 ^- + e<f) £ w on both sides of 



(I3.37p . we obtain 

/ \V6 £ \ 2 dx < A 2 / \eV(N j 
Jn * Jn 



dxi 



■F)\* dx + £) 2 I \eV(w6 £ )\ 2 dx 



A 



(3.46) 



where 4> £ is a cut-off function, satisfying 



(jf = 1, on <9f2 
< (f) £ < 1, on n £ 
\eV(j) £ \ < C 

£ e c°°(n), 



(3.47) 



with Vt £ = {x eVt\ dist(x, dVL) < e}. 

For Ii, we directly use the result in [21 



du°, 



/i<^(l+||iV||^ (y) ) 2 || — ||^ 1(n) 



(3.48) 



Please note that the above estimate only depends on the ||iV(y) rather than 

11-^(2/) || wi.°°(y) • This is the contribution of Suslina 
For I2, we have 



J 2 < e 2 \V x w(x, y)\ 2 dx + / \V y w{x 1 y)\ 2 dx + / \w(x, y)\ 2 dx 
Jn £ Jn E Jn e 

< Ce||'u;|| 2 yi, 00 ( nxY ). 
By ([Sag), (EHBj) . (|328D, and Q333D, we complete the proof. 



(3.49) 



□ . 



Remark 3.3 Please note that the estimate in l[3.45\ ) is of order e. The final estimate 
( 13 . 1 3[) decays to the order , due to the oscillation of the corrector term 9 s ( 13 . 3 7|) at 
the boundary. For the one-dimension problems, the boundary decays to isolated points 
and this kind of oscillation at the boundary does not appear any more. The estimate 
( 13 . 13[) can be improved to 



(u £ 11^(0,!) <Cfe- (3.50) 
Remark 3.4 As to the convergence of potential energy, we have as e — > 

| / (Vu £ • A £ Vu £ + F £ ■ Vu £ ) dx - [ (Vu° • A°Vu° + F° ■ Vu°) dx\ 0. (3.51) 
Jn Jn 

So far we have establish the mathematical homogenization theory for (I3.4j) . we now 
come back to the physics problem ( 13. ip . 

By Theorem 13 . 2 1 and unit transformation, we get the homogenized problem of ( 13. II) . 



(3.52) 



(3.53) 



-V ■ {K°V Po {x)) = f°(x) + V ■ F°(x), in Q D = (0, D) d , 
Po = 0, on dVlo, 

where 

p (x) = D 2 p (x) 
K° = K° 

with K°, f°, and F° defined as follows 

«5 = (^(»)+^©y, 

f°(x) = (f(x,y)) Y , (3.54) 
F°(x) = (F(x,y) + K(y)V y w(y)) Y , 

and Ni(j) = N^(y), w(x, f) = w(x,y) solving the cell problems 

-V, ■ (K(y)V y W{y)) = V, ■ (K(y)e>), in Y, 

N j (y) is periodic in Y, (N^) Y = 

-V, • (K{y)V y w(x,y)) = V y ■ {F(x,y)), in Y, 

w(x,y) is periodic in Y, (w)y = 0. 
By unit transformation, (13.531) . and (13.131) . we can obtain the next theorem 



(3.55) 



(3.56) 



Theorem 3.5 If p is the solution of ( 13. ip . p\ is defined as follows 

Pl =p + lN^ + lw, (3.57) 

dxj 

with po solving (I3.52j) and f(x,y), F(x,y) are bounded, smooth enough and Y-period 
with respect to y, then there exists a positive number C independent of D, such that 

^ L |v( £ - § )|2<fa s 4- (3 - 58 » 



D 



As to the convergence of potential energy, we have 

\E(p) - E (po)\ ^ asD oo, (3.60) 

*w-^£(vg)*v(S) + £v(S)) fc 

Above theorem explains in what sense we can accept the homogenized problem ( 13.521) . 

Remark 3.6 In one-dimension case, by Remark \3.3[ our results are changed to be: 

£f <-|)l^< C (i) 2 , (3.61) 



with 



\E(p) - E (p ) \ ^ as D oo, (3.63) 



4 Example 



In this part, we show a one-dimensional example to verify our results . We consider 
the following problem: 

-d x (Kd xP ) = -l + d x ((x + % + C)K) m(0,D), 

(4.1) 

P(0) = p(D) = 0, 

, ^ 1 / . (2ixD\ ( I V 1 ^2tt/A (I 

where K = -s — - and C = — sin — - — + — — cos 



2 + cos(^) " 2tt \ I J \2tt J D \ I J \2n J D 

It is clear that the source term has micro-structure, which we have discussed in §3. 
The solution is 

x 2 D I . f2nx\ ( I \ 2 (2-nx 



p = x + x — sin — — I + — cos 



2 2 2tt V I J \2n J V I 

I ^ 2 



-Cx-[-) . (4.2) 



The homogenized problem has the form as 

-d x (K Q d xPo ) = -l + d x (f ) , in (0, D), 

Po (0) = Po (D) = 0, 

x 2 D x _ sin(2f ) 



(4.3) 



(4.4) 



with K° — \ and the solution p = — — —x. N(-) = — — - — is the solution of cell 
problem 

-d y (Kd y (N + y))=0, inF = (0,l), 
N(y) is Y-periodic with respect to y and (N(y)) Y = 0, 

x d s i n (2zE£) 

and w(x, — ) = (x H h C) — — - — is the solution of cell problem 

I 2 An 

-d y (K(y)d y w(x,y)) = d y (F(x,y)), in Y = (0, 1), 
w(x,y) is Y-periodic with respect to y and {w(x,y))y = 0. 

By direct computation, we obtain 

/ (p-Pofdx 
Jo 

D/ i f^A f l Y f 2 ™\\\ f D ^,2, f D ( i i - 

x— sin I —r- ) + ( — ) cos I —r~ ) | ax + (Cx) ax + ( — ) ax 

Jo Jo 



(4.5) 



o 



l3 



2vr 1 V 2tt 7 \ I .In .In V2tt 



=D a — — + —) — D \ — — sin — — + D ( ^-^ - cos [ 



4tt£>\ Z 2 



8tt 2 



/ Z 5 /4tt£)\ / 3 , (A-kDX. 
+ (T28^ Sm (^J"^ Sm l^J, ) - 



So we have 



D 

which is consistent with the theoretical result (I3.62p . 
If we denote by 



then we have 



Pi = Po + lNd x p + Iw 

I . f2ux\ IC . f2nx\ 
= Po + x-sm I— +-sm — ) , (4. 



I V {2irx\ ( I \ 2 IC (2-kx 



p- Pl = (- cos — - - -Cx + -sin( — ). (4.9) 



D 2 . f° f I . f2nx\ „ C / / 2ttx x x 2 

JO 



(5 x (p-pi)) cte = / ( — — sin I — ^— J — C + — cos f — ^— ) ) dx 



^ + ^MS-3^> 8i »l^).(4.10) 



(A<kD\ 



We obtain 

This is also consistent with the theoretical result (I3.6ip 
As to the convergence of the energy, we have 

D pD 

d x p ■ Kd x pdx + / Fd x pdx 
(I Jo 
D 

fpdx 







D 3 CD 2 / Z 2 /2vrD\ I 2 \ I . (2nD\ u 1 . 

= h £> — o cos H sin — — , (4.12) 

/ a xPo • K°d xPo dx + / F°d xPo dx = — , (4.13) 

Jo Jo 

so we have 

lnM „ . C 1 / / 2 /2tt£>\ Z 2 \ l (2irD\ 

\E(P) - E M \ = - + - I — cos I — \ + — \ - — sm ( — ) - 



as D — > oo , with 



E (po) = 




This is also consistent with the theoretical result (13.631) . 
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